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Ram pressure stripping in a galaxy formation model. I. A novel 
numerical approach 
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ABSTRACT 

We develop a new numerical approach to describe the action of ram pressure stripping (RPS) 
within a semi-analytic model of galaxy formation and evolution which works in combination 
with non-radiative hydrodynamical simulations of galaxy clusters. The new feature in our 
method is the use of the gas particles to obtain the kinematical and thermodynamical proper- 
ties of the intragroup and intracluster medium (ICM). This allows a self-consistent estimation 
of the RPS experienced by satellite galaxies. We find that the ram pressure in the central 
regions of clusters increases approximately one order of magnitude between z = 1 and 0, 
consistent with the increase in the density of the ICM. The mean ram pressure experienced 
by galaxies within the virial radius increases with decreasing redshift. In clusters with virial 
masses M v ; r ~ 10 15 hT 1 M©, over 50 per cent of satellite galaxies have experienced ram pres- 
sures ~ 10~ n h 2 dyn cm~ 2 or higher for z < 0.5. In smaller clusters (M v j r ~ 10 14 hr 1 M ) 
the mean ram pressures are approximately one order of magnitude lower at all redshifts. RPS 
has a strong effect on the cold gas content of galaxies for all cluster masses. At z = 0, over 
70 per cent of satellite galaxies within the virial radius are completely depleted of cold gas. 
For the more massive clusters the fraction of depleted galaxies is already established at z ~ 1, 
whereas for the smaller clusters this fraction increases appreciably between z = 1 and 0. This 
indicates that the rate at which the cold gas is stripped depends on the virial mass of the host 
cluster. Compared to our new approach, the use of an analytic profile to describe the ICM re- 
sults in an overestimation of the ram pressure larger than 50 per cent for z > 0.5. 

Key words: galaxies: clusters: general - galaxies: clusters: intracluster medium - galaxies: 
evolution - galaxies: formation - intergalactic medium. 



1 INTRODUCTION 

Comparison between galaxies in clusters and in the field reveals 
remarkable differences in their physical properties, which suggest 
the important influence of environment on galaxy formation and 
evolution. Large galaxy surveys have shown that th e colour distri- 
butio n of the galaxy pop ulation is clearly bimodal dStrateva et al.l 
l200ll ; iBaldrv et al]|2004l) . One peak is formed by red, non star- 
forming galaxies, while the other is populated by blue galaxies with 
active star formation. In the local Universe, the fraction of red se- 
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quence galaxies at a given stellar mas s depends on environment 
dBalogh et a l. 2004; B aldrv et alj|2006f) . with galaxy clusters hav- 
ing larger fractions of red, early-type galaxies than the field. Envi- 
ronment also influences the star formation rates (SFRs) of galaxies, 
as these are strongly suppressed in dens e environments for galax ies 
over a large range of stellar masses (e.g. lKauffmann et al.|[2 004). 

Another quantity which shows dependence on environment is 
galaxy morphology. Spiral galaxies tend to be r arer in the cen- 
tral regions of clusters than ear ly-type galaxies toresslerl fl980; 
IWhitmore. Gilmore. & Jones 1993), and cluster spiral galaxies dif- 
fer from those in the field in several characteristics that can 
be correlated with environment (for an extensive review see 
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IBoselli & Gavazzill2006l) . In clusters, spiral galaxies are not only 
redder but also more HI deficient than similar galaxies in the field, 
with this deficiency increasing for galaxies closer to the clus- 
ter centre (Havnes, Giovanelli, & Chincarinil [l98i | Solan es et al .1 
1200 ll ; iHughes & Cortesell2009l) . In field spirals the gaseous discs 
typically extend beyond the optical radius, but the opposite trend 
is seen in clusters, where a significant fraction of spiral s have 
truncated gaseous discs (e. g. Koopm ann & KennevI 120041 1 20061 ; 



iKoopmann. Havnes. & Catinellall2006l) 

A likely scenario is that galaxies are transformed from blue 
star-forming systems into the red and passive population, but it is 
still unclear which physical processes are the key ones responsible 
for this transformation, and what their relative importance may be. 
The observed differences between cluster and field galaxies out- 
lined above suggest that the transformation could be the result of 
processes that remove gas from galaxies, suppressing their star for- 
mation, and/or alter the galaxy morphology, both effects happening 
preferentially in high-density environments. Galaxy-galaxy inter- 
actions and mergers are a mong these possible mechanisms as su g- 
gested by some theoretical ( Moore et al .1 19961 : Pere z et al .l2006allbh 
and obs ervational (e.g.lEllison et al. 2008, 201o| ; |Perez et alj2009h 
results. iMoore et al. ( 1996, 1999) proposed that the evolution of 
galaxies in groups and clusters is governed by the combined ac- 
tion of multiple close encounters between galaxies and the tidal 
interaction with the group potential, a process that they call 'galaxy 
harassment'. In clusters, however, the relative velocities of galaxies 
are higher than in galaxy group s and, consequently, the interaction 
times are shorter jMihosll2004f) , so this process is expected to be 
more significant in less dense environments. 

Another process that might be relevant is the removal of 
the hot diffuse gas halo of a galaxy after its infall into a group 
or cluster. This has been called 'st rangulation' or 'starvation' 



dLarson. Tinslev. & Caldwell! Il980l : iBalogh, Navarro. & Morris] 
l200d : iKawata & MulchaevT 2008). The removed gas becomes 
part of the overall intragroup or intracluster medium (ICM), the 
affected galaxy cannot accrete any more gas via cooling flows, and 
with a moderate SFR it will consume all of its cold gas within a 
few Gyr, ending its star formation and becoming gradually redder 
as its stellar population ages. Strangulation is a standard ingre- 
dient in most semi-analytic models of galaxy formation, which 
successfully reproduce observed global properties of galaxies such 
as luminosity functions, colour distribution s and mass-metallicity 
relations (e. g . iBaugh. Cole, & Frenkl 1 19961: ISpringel et alj I2OO1I; 
Bower et alj 120061: ICroton et alj 12004 iLagos, Cora. & Padillal 



2008L hereafter LCP08). As commonly implemented, as soon as a 



galaxy becomes a satellite, its hot gas halo is assumed to be shock 
heated to the virial temperature of the group and immediately 
removed from the galaxy. This effect proceeds in the same fashion 
in groups of all masses. 

Ram pressure stripping (RPS) of the cold gas in galac- 
tic discs could also play an important role. Galaxies in clus- 
ters move through the hot diffuse ICM gas at velocities 
that could be close to supersonic l lFaltenbacher & Diemana 
20061). and s o will experience considerable ram pressure (RP). 



Gunn & Got3 d 1971 hereafter GG72) proposed that when the 
RP exceeds the gravitational restoring force of the galaxy, 
its cold gas will be pushed out. In recent years, RPS 
has been extensively studied using hydrodynamical simula- 

owen | 19991; 
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lias uccn extensively stuuicu using iiyuiuuyiiaiincai ^ 
tions of individual galaxies (e .g. Abadi, Moore, & Bower 
Ouilis, Moore. & Bowed |2000|; iMarcolini. Brighenti. & D' 



20031 : iRoediger & Bruggen 20061. l200l iKronberger et alj |2008|) 



which suggest that the GG72 estimate is a fairly good approx- 



imation in most situations, and that the time-scale for gas re- 
moval is ~ 10 - 100 Myr. RPS acts only on the gaseous com- 
ponents of the galaxy, so its characteristic signature is the pres- 
ence of truncated gas discs while the stellar discs remain un- 
altered. In deed, there are se v eral observat i onal c andidates for 
RPS ( e.g. ICrowl et alj 120051 : IBoselli et alj 120061: [Cortese et al.l 
l2007l : ISun! Donahue. & Voill l200l I Vollmerl 120091) . RPS also 
seems to be relevant for dwarf ga laxies in less den s e en- 
vironments, such as galaxy groups dMcConnachie et al. 1 120071: 
IMastro pietro. Burkert, & Moore 2008). A recent revie w of both 
simula tions and observations of RPS can be found in IRoediger 
d2009h . 

All the physical processes described above are not mutually 
exclusive. In a hierarchical universe, many galaxies in clusters were 
previously members of smaller systems where a combination of 
strangulation, harassment and/or RPS could start to suppress star 
formation ( t his has been called 'galaxy pre - processing'; see e.g. 
lFuiitall2004lMihosll2004l ; ICortese et all2006l : |Perez et alj|20091) . 

Semi-analytic models of galaxy formation include galaxy 
mergers and strangulation as standard elements, but the effect 
of RPS has been included i n such models only in few cases 
(Okamoto & Nagashima 2003; Lanzoni et al. 2005) where no sig- 
nificant influence of RPS was found on the analysed galaxy proper- 
ties. These previous works use dark matter (DM)-only simulations 
to generate the merging history trees of DM haloes, which are then 
used by the semi-analytic model to generate the galaxy population, 
and the properties of the ICM and velocity distributions of galaxies 
are modeled using analytical approximations. 

A more recent work by IBriiggen & De Lucial (2008J, here- 
after BDL08) studies the distribution and history of the RP expe- 
rienced by galaxies in clusters, comb ining a semi-ana lytic model 
with the Millennium Simulation dSpringelet ail 12005): they use 
the properties of the DM particles to track the positions and ve- 
locities of simulated galaxies. Again, since the Millennium is 
a DM-only simulation, BDL08 have to resort to analytic mod- 
els for the ICM properties, assuming the hypothesis of hydro- 
static equilibrium for the gas w ithin a DM halo described by the 
iNavarro. Frenk. & White! 1 19971 hereafter NFW) profile. The lat- 
ter is a good approximation to the DM profiles if the haloes are in 
dynamical equilibrium; this is likely not the case for high-redshift 
haloes , where the bulk o f the star formation activity takes place 
(e.g lMadauetal]|l996l : Iciardi. Stoehr. & White! [20031 : iHopkinsI 
2007). Since the dynamics of the I CM may play an important role 
1 Sun vaev. Norman. & Brvanll2003l) and certainly the hypothesis of 
hydrostatic equilibrium may not hold for all haloes, our aim is to 
develop an improved model for RPS which takes into account the 
ICM dynamics through simulations that include gas physics. 

We model RPS by adopting the criterion proposed by GG72 
and by implementing this process in the semi-analytic model 
of galaxy formation and evolution Semi-Analytic Galaxies (sag; 
LCP08), which is combined wi th hydrodynamic al cosmological 
simulations of galaxy clusters l lDolag et al JI20051) . The novel fea- 
ture of our implementation is the fact that the thermodynamical 
and kinematical properties of the ICM, which are involved in the 
estimations of the RP experienced by each galaxy, are provided by 
the underlying simulations. In this paper, the first of a series, we 
focus on the study of the distributions of RP experienced by satel- 
lite galaxies in clusters of different masses, the dependence of RP 
with clustercentric distance, and how these quantities evolve with 
redshift. We compare our results with those obtained by assuming 
NFW profiles for the density distribution of the ICM. A forthcom- 
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ing paper will deal with the influence of RPS on galaxy properties 
such as luminosities, colours and star formation histories. 

This paper is organised as follows. In Section [2] we briefly 
describe the semi-analytic model sag used in this work and give 
details of the simulations used. Section [3] contains a detailed de- 
scription of the way in which RP is estimated from the hydrody- 
namical simulations, and the implementation of the RPS process in 
the semi-analytic model. Section|4]presents the distributions of RP 
that we obtain from our cluster simulations and compare them to 
the results obtained by adopting a NFW profile for the ICM under 
the hypothesis of hydrostatic equilibrium. We then analyse the ef- 
fects of RPS on the gas content of galaxies. Finally, in Section[5]we 
summarize our main findings. 



2 NUMERICAL SIMULATIONS AND THE 
SEMI-ANALYTIC MODEL 

To develop our RPS model we use a hybrid numerical approach 
which combines cosmological non-radiative TV-body/smoothed 
particle hydrodynamics (SPH) simulations of galaxy clusters with a 
semi-analytic model of galaxy formation and evolution. In the cos- 
mological simulation, the DM haloes and their substructures are 
identified and followed in time, building up detailed merger trees 
which are used by the semi-analytic code to generate the galaxy 
population. The main advantage of this hybrid procedure is that the 
use of semi-analytic codes allows the exploration of a larger dy- 
namical range than fully self-consistent simulations, at a fraction 
of the computational cost. 

In the present work we implement the RPS process in the 
semi-analytic code developed by LCP08. We take the kinematical 
and thermodynamical properties of the hot diffuse gas which per- 
meates DM haloes directly from the underlying SPH simulation, 
thus avoiding the use of analytical approximations. In this section, 
we briefly describe the simulations and the semi-analytic model on 
to which we graft our RPS scheme described in detail in Section[3] 



2.1 Non-radiative iV-body/SPH simulations of galaxy 
clusters 

In this study, we use the high-resolutio n hydrodynamic non- 
radiative simulations of galaxy clusters by |Polag"e t al. (2005]). The 
simulated clusters were originally extracted from a DM-only simu- 
lation with a box size of 479 h^ 1 Mpc of a flat ACDM model with 
n m = 0.3, Qa = 0.7, H = 100 /iT 1 km s _1 Mpc -1 where h = 0.7, 
a baryon densit y Slj, = 0.039 and a normalization o f the power spec- 
trum as = 0.9 dYoshida. Sheth. & Diaferid l2001). The Lagrangian 
regions surrounding the selected clusters have been resimulated at 
higher mass and force resolution using the 'zoom ed initial condi- 
tions' technique ( iTormen. Bouchet. & White! 19971) . Gas was intro- 
duced in the high-resolution region by splitting each parent particle 
into a gas and a DM particle. The mass resolution is the same for all 
simulations, being 771dm = 1.13 x 10 9 h -1 M@ for DM particles 
and m g!ls = 1.69 x 10 s ft -1 Mq for gas particles. However, the 
identification of DM haloes was based only on the DM particles, 
with their mass increased to its original value. 

The simulations have been carried out using gadget-2, a par- 
allel Tree-SPH code with fully adaptiv e time-step ping and explicit 
conservation of energy and entropy (Springel 2005). For the force 
resolution, the gravitational softening is fixed at e = 5 h~ kpc 
in physical units at z ^5, and for higher redshifts it shifts to 



e = 30 h~ kpc in comoving units. For each simulation, 92 snap- 
shots were stored between redshifts z = 60 and z = 0. The simu- 
lations considered here include only non-radiative physics and the 
original formulation of artificial viscosity within SPH; this is a rea- 
sonable approximation for simulations of clusters, since the bulk of 
the halo gas has very long cooling times. 

From these simulations, we have selected eight 
cluster-sized haloes divided in two sets: five clusters with 
Mvk — 10 14 H~ 1 Mq (hereafter, G14 clusters) and three clusters 
with M vir ~ 10 15 /i _1 M Q (G15 clusters). These DM haloes are 
first identified in the simul ation outputs by means of a friends- 
of-friends (fof) al gorithm dDavis et al 1 1985b . Subsequently, the 
subfind algorithm JSpringel et alj 1200 II) is applied to the haloes 
detected by fof in order to find self-bound DM substructures which 
we call subhaloes. To build the trees, we extract from the simula- 
tions all subhaloes consisting of 10 or more D M particles, since 
small er ones are usually dynamically unstable ( Kauffmann et al. 
1999). Additionally, SPH particles in the simulations provide the 
density, spatial and velocity distributions of the intergalactic and 
intracluster media. 

2.2 Semi-analytic model of galaxy formation 

The DM halo merger trees and the information from the SPH par- 
ticles are then used as input for the semi-analytic mod el of galax y 
formation sag (LCP08), based on the one described bv lCorai ( 2006), 
which includes the effects of radiative cooling of hot gas, star for- 
mation, feedback from supernovae explosions, chemical enrich- 
ment and galaxy mergers, sag includes the effect of strangulation 
as described in the Introduction. LCP08 updated it to include black 
hole (BH) growth and feedback from active galactic nuclei (AGN). 
In addition to this major improvement, sag allows starbursts in three 
different ways: both in major and minor m ergers, and w hen disk in- 
stabilities occur. The reader is referred to lCoral (2006) and LCP08 
for the full details of these implementations. In the present work, 
this model is further modified to include the effect of RPS on galax- 
ies due to the hot intergalactic gas as described in the next section. 

The galaxy catalogue is built up by applying the semi-analytic 
model to the detailed DM subhalo merger trees extracted from the 
hydrodynamical simulations. Similarly to other semi-analytic mod- 
els, in this 'subhalo scheme' arising from the identification of DM 
substructures within fof groups, the largest subhalo in a fof group 
is assumed to host the central galaxy of the group, located at the 
position of the most bound particle of the subhalo. These galaxies 
are designated as central or type galaxies, and each fof halo has 
only one. Central galaxies of smaller subhaloes contained within 
the same fof group are referred to as halo or type 1 galaxies. The 
subhaloes of these galaxies are still intact after falling into larger 
structures. There is a third group of galaxies generated when two 
subhaloes merge and the galaxy of the smaller one becomes a satel- 
lite of the remnant subhalo: these galaxies are called type 2 galax- 
ies. A type galaxy can have satellites of types 1 and 2 and, fur- 
thermore, a type 1 galaxy may itself have type 2 satellites. In the 
following, whenever we use the term 'satellite galaxies' we will 
be referring to both type 1 and 2 galaxies. We assume that type 2 
galaxies merge with their corresponding subhalo central galaxy on 
a dynamical friction time-scale. 

For each of these galaxies, sag provides information on the 
stellar mass, cold disc gas, hot gas within DM haloes, BH mass, 
AGN activity, star formation histories, magnitude s in several bands 
accord ing to the stellar population models by iBruzual & Chariot! 
J2003I) and chemical abundances of the different baryonic compo- 
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nents. Galaxy positions and velocities are traced by tracking the po- 
sition and velocity of the most bound particle of their host DM sub- 
halo; for type 2 galaxies, we use the most bound particle identi- 
fied at the last time there was a substructure (as in BDL08). In a 
forthcoming paper we will test how the effect of RPS depends on 
the choice of alternative recipes for the different physical processes 
such as the star formation scheme and the galaxy orbits. 



3 RAM PRESSURE STRIPPING OF COLD DISC GAS 

In this Section we describe our self-consistent implementation of 
RPS in the semi-analytic co de. Our RPS model is based on the 
simple criterion proposed bv lGunn & Gottl dl972[) : the cold gas of 
the galactic disc located beyond galactocentric radius R will be 
stripped away if the RP exerted by the ambient medium on the 
galaxy exceeds the restoring force per unit area due to the gravity 
of the disc: 



Picmv 2 5s 27rGE diS c(i?) Ecoid (R). 



(1) 



Here picM is the ambient density of the ICM at the current posi- 
tion of the galaxy, v the velocity of the galaxy relative to the ICM 
and Edisc, Scold are the surface densities of the galactic disc (stars 
plus cold gas) and of the cold gas disc, respectively. Although the 
description of RP phenomena is usually quoted in the context of 
galaxy clusters, it is valid in general, and we apply this criterion to 
the hot gas contained within all DM haloes at all redshifts. 

In order to determine the RPS, we need an estimation of each 
of the parameters involved in condition |QJ, namely the properties 
of the ICM on the left-hand side and the scalelength of the galactic 
discs necessary to evaluate the right-hand side. Below, we explain 
the strategies used to estimate these quantities within our hybrid 
model. 



3.1 Determination of the ICM properties 

As already mentioned in the Introduction, the main difference be- 
tween ou r RPS implementation in a semi-analytic model and previ- 
ous ones dOkamoto & Nagashima 2003; La nzoni et al.l2 005) is that 
we use the fhermodynamical properties of the ICM derived from 
the underlying hydrodynamical simulation. For all satellite galax- 
ies in the simulations, we determine the local ICM density picm and 
velocity v relative to the ICM by searching for all gas particles in 
a sphere of radius r ng b oc r v i r centred on each galaxy, where r v i r is 
the virial radius of the host DM subhalo of the satellite. For type 2 
galaxies, the value of r v j r corresponds to the last time the galaxy 
was identified as either type or type 1 . If less than 32 particles are 
found, we take the closest 32 neighbours instead. 

Some of the particles selected by this procedure might be grav- 
itationally bound to the DM subhalo. In such a case, ICM densities 
would be overestimated, and the velocity of the galaxy relative to 
the mean ICM would artificially be v ~ 0. Thus, these particles 
need to be discarded since we are interested in estimating the prop- 
erties of the ambient medium which is responsible for the RPS. To 
do this, once all gas neighbour particles have been found, the higher 
density ones are filtered out by using an iterative procedure: we de- 
termine the median gas density p m of the selected particles, discard 
all particles with p > / lm p m and repeat this procedure until the 
median density converges. We find that choosing r ng b = 2.5 r V k and 
frm = 2, a smooth ICM density distribution can be recovered, re- 
moving substructure in the gas without significantly affecting the 
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Figure 1. Top: mean ICM density profiles of the simulated clusters at z = 0. 
Open diamonds and filled squares show values for the G14 and G15 clus- 
ters, respectively. Data points for the G15 clusters are shifted to the right 
for clarity. Dashed lin es show the average density profiles determined by 
IVikhlinin etal] <2006l) for a sample of nearby relaxed clusters. Bottom: re- 
lation between disc scalelength and /-band luminosity for disc galaxies, 
for a ran of the SAG model that includes the RPS effect. We selected galax- 
ies with Afbaryon 10 8 h~ x Mq, Mb u i gc /M s , e ii ar ^ 0.95 and with central 
surface brightness po(B) < 22 mag arcsec" 2 . For this plot we take all the 
galaxies in one of the largest simulation boxes, including both cluster-sized 
haloes and smaller systems. Solid and dashed lines show the mean observed 
relation and la scatter from Courteau et al. 1 20(53). 



median density profile. The final step is then to determine the veloc- 
ity v of the galaxy relative to the mean motion of the gas particles 
that remain after filtering. 

By calculating the ICM properties in this fashion, we obtain a 
self-consistent method which does not introduce any additional free 
parameter into the model and does not need to assume that the ICM 
gas is in hydrostatic equilibrium. In fact, it automatically takes into 
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account any local variation of the density and/or the velocity field 
due to the dynamics of the gas. The method is not overly sensitive to 
the values of f m and r ng b chosen; r ng b is a compromise between the 
need to consider a radius large enough to obtain a fair sample of the 
local environment, but not so large that the search spans a large por- 
tion of the cluster volume. With / rm > 3 the method becomes inef- 
fective at removing substructure, and choosing / rm ^ 1 the method 
underestimates the median ICM density when compared to obser- 
vations. 

The density profiles we obtain from this procedure, for the 
case of the cluster-sized haloes, are in excellent agreement with 
those deter mined from X-ray observations of ga l axy clusters at 
z ~0 (e.g. lSchindler. Binggeli. & B6hringeilll999l ;l Vikhlin in et al.1 
2006). This can be seen in the top panel of Fig.[T] where we com- 
pare the mean ICM density profiles obtained for the sim ulated 
clusters with the profiles obtained by I Vikhlinin etai] ( 120061) for a 
sample of nearby relaxed clusters, which span a temperature range 
0.7 — 9 keV that comp rises the mass-we ighted temperatures of the 
simulated clusters (see lDolag et alj2005l) . 

3.2 Sizes of galactic discs 

In the semi-analytic code sag, the gas acquired by a galaxy via cool- 
ing flows is assumed to settle onto a thin exponential disc with sur- 
face density E(ii) = Eo exp(— R/Rd), where Eo is the central 
surface density and Rd the disc scalelength. If the outer disc ra- 
dius Rase S> Rd, then these quantities are related to the total disc 
mass Md through 



2wR 2 / 



(2) 



To determine the disc s cale length we follow the model developed 
bv lMo. Mao, & White! d 19981 , hereafter MMW). Assuming a NFW 
profile to represent the DM distribution, the MMW model allows 
the calculation of Rd given the halo virial mass M y - a , the mass 
fraction of baryons that settle on to the disc md = Md/M v t, the 
halo concentration c and the dimensionless halo spin parameter 



A = J|£| 1/2 G _1 M 



-5/2 



where J and E are the angular momen- 



tum and energy of the halo, respectively. 

The MMW model assumes that the DM halo responds adia- 
batically t o the slow assembly of t he disc, following the standard 
scheme of Blumenthal et al.1 dl986h . The initial NFW mass profile 
M(ri) and the profile after contraction Mf(r) are then related by 



M f (r) = M d (r) + M b + M(n)(l - m d - m b ), 



(3) 



where Mt, is the total bulge mass (considering bulges as point 
masses, for simplicity), m b = M b /M v i, , and M<j(r) is the disc mass 
within r given by 



M d (r) = 2-kT.qRI 



1 



Rd 



exp(-r/7? d ; 



(4) 



The implicit assumption is that the initial distribution of baryons 
parallels that of the DM, and those that do not end up in the disc or 
in the bulge remain distributed in the same way as the DM. These 
considerations allow MMW to express the disc scalelength as 



Rd = -7= ( — J \rwfc 1/2 /fl(A,c,m d , j d ), 
V2 \m d J 



(5) 



where jd = J d /J, with Jd the total disc angular momentum, and 
fc fa are factors given by 

2 / c \0.7 

3 + V2L5 



fc 



(6) 



f R (X,c,m d ,j d ) = 2 



, 2 V c (R d u) ' 
u — — du 



V'vii 



(7) 



(see MMW for the full details). If the specific angular momentum 
of the disc material is the same as that of the parent halo, the n 
the ratio jd/rrid should be close to unity dFall & E fstathiou 19801) . 
Thus, we set jd = md and assume the angular momentum of the 
bulge to be negligible. 

For a given set of values of V y \ r , c, A, m d and jd, equations 
l[3} - (0 must be solved by iteration to yield Rd - However, including 
this iterative method in the semi-analytic model is extremely time 
consuming. Most semi-analytic models use a fitting formula for Rd 
determined by MMW (their equation 32), but this is valid only for 
the case of pure disc galaxies and with an accuracy of 15 per cent 
quoted only for values of md in the range 0.02 < mj < 0.2. We 
find that galaxies in our simulations cover a much broader range 
in md', besides, we also want to take into account the presence of 
bulges and their gravitational effect on the size of the final disc. 
Hence, we use the iterative procedure to generate a set of look-up 
tables for each snapshot of the underlying simulation and for a grid 
of values of md, m b , c and A. These tables are then used by the 
semi-analytic model to find Rd by interpolation. 

We assume that the process of disc formation takes place only 
for central galaxies. When a galaxy becomes a satellite, the accre- 
tion of cold gas is halted due to strangulation, so its disc cannot 
grow any further. We then take Rd to be frozen at the value it had 
at the last time the galaxy was identified as a central. Additionally, 
if a galaxy suffers a major merger or if the galactic disc becomes 
unstable, the remnant is a spheroid and so Rd is set equal to zero 
after such an event. 

The values of mj and m b are calculated for each central 
galaxy from the properties given by the semi-analytic code. The 
concentration of their host DM haloes are determin ed as a func- 
tion o f M v i r and redshift by following the model bvlBullock et al] 
d200lh . with the parameters quoted bv lWechsler et alj J2006t) . The 
spin parameter A is determined for each halo by computing its 
total energy and angular momentum. Studies based in iV-body 
simulations have found that the distribution of A for DM haloes 
can be well described by a lognormal functio n with mean 0.03 < 
Ao 5, 0-05 and dispersion a ~ 0.5 (e.g. IWarren et alj 1 19921 : 
iBett et alfcOOllMaccio et alj2007l) . For the simulations used here 
Ao = 0.042, but the distribution deviates from the lognormal for 
high values of A. This happens because the mass resolution in N- 
body simulations affects the determination of angular momenta, in- 
dependently of the method used to find DM structure, ge nerating an 
artificial high-A tail in the distribution < lBettetal.tt2 007'l. The value 
of the spin parameter correlates with the degree of equilibrium of 
the DM halo, with haloes which are far from being virialized hav- 
ing large val ues of A. As a con sequence, we use the simple criterion 
proposed bv lBett et ai]( l2007h to filter out these anomalous values. 
For each halo, we calculate an 'instantaneous virial ratio' 



Q= \2T/U + l\ 



(8) 



and consider that haloes with values Q < 1 are in a quasi- 
equilibrium state. This allows us to correct the distribution of A for 
mass resolution effects. For central galaxies of anomalous haloes 
(Q > 1), we set R d = A r vn /V2 (MMW, their equation 12). 

This procedure generates a distribution of sizes of galactic 
discs which is in good agreement with observations. The bot- 
tom panel of Fig. Q] shows the relation between disc scalelength 
and 7-band luminosi ty, compared to the observed relation from 
ICourteau et alj (|2007). The data were extracted from a run of the 
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sag model which includes the RPS effect as described in Sec- 
tions l3.1l and l3,3l (see also Section l4~2l l. The plot includes all galax- 
ies with A/baryon 10 8 Mq within one of the largest sim- 
ulation boxes, comprising not only cluster-sized haloes but also 
smaller systems. We select as spirals those galaxies with central 
surface brightness fio(B) < 22 mag arcsec -2 to e xclude low 
surface brightness (LSB) systems dlmpev et al.l 1200 lT) . and with 
r, h , cs h = Mbuig e /M ste nar < 0.95 (following LCP08). Varying r, hreS h 
between 0.5 and 0.95 does not result in any significant change in 
the distribution of morphological types with galaxy stellar mass. 

The agreement between the model and observations is very 
good both in terms of the slope and the scatter of the relation. 
For a given luminosity, galaxies with higher values of A (i.e. with 
higher angular momentum) have larger discs. Once the artificially 
large values of A have been corrected, very few galaxies (< 5 
per cent) fall in the region above the upper dashed line in Fig. Q] 
and those are mostly LSB galaxies. Disc galaxies are not found in 
the region below the lower dashed line (and below the mean for 
Li > 10 11 h~ 2 LQ) because such discs become unstable. 



3.3 Modelling ram pressure stripping 

With the relevant properties of the ICM estimated as previously 
described we can now construct our RPS model. The mass loss due 
to RPS is calculated analytically as follows. For exponential discs 
of stars and gas, it can be shown from condition Q} that RPS will 
remove from a galaxy all the gas beyond a stripping radius i? st i 
given by 



-0.5R d In 



PICM V 



2ttGE 

,disc£o,c 



(9) 



Here Eo.coid and Eo.disc are the central surface densities of the cold 
gas disc and of the galactic disc. We assume, for simplicity, that 
both disc components (stars and cold gas) have the same scale- 
length. To avoid extremely large values of _R stl , which may lead 
to divergences in the calculations, we ignore the effect of RPS if 
Rsa/Rd ^ 20; in such cases the fraction of gas removed is negligi- 
ble. 

Our calculation of mas s loss follows that of iLanzoni et al.l 

d2005l) . The previous study bv lOkamoto & Nagashimal J2003h used 
a much rougher estimation of RPS: for each galaxy they calcu- 
lated a mean gravitational restoring force, and if the RP exceeded 
its value, all the cold gas was removed; if the RP was smaller, 
no gas was lost. Thi s is a crude assumption, and we consider the 
ILanzoni et al.l d2005h approach a better representation of the contin- 
uous gas removal caused by RPS. A difference between our work 
and ILanzoni et alj d2005h is that they include a factor cos 2 i in the 
left-hand side of (TJ to account for the disc inclination, where i 
is the angle between the normal to the plane of the disc and the 
direction of motion of the galaxy through the ICM. We have cho- 
sen to neglect this factor because studies based on SPH simula- 
tions agree that the galaxy inclination does not have a significant 
influence on the amount of removed gas; strong RPs that strip a 
face-on galaxy also strip a galax y moving edge- on, although on 
a slightly longer time-scale (see iRoedigen I2009L and references 
therein). We have also chosen to ignore the possibility of en- 
hanced star formation as a result of RP compres s ing the molecu- 
lar ga s of the galaxy I Fuiita & Nagashimal 19991; Bekk i & Couch] 



120031; iKronberger et alj 120081 ; iKapferer et alj 120091), since ob 
servational evidence is contradictory dKoopmann et al.l 1 20061 ; 
iKoopmann & KenneV 2006; Abramson & Kenney| |2009l) and the- 



oretical work needs confirmation due to the complex nature of star 
formation. 

In sag the evolution of the galaxy population is followed by 
solving differential equations describing the physical processes 
involved, such as gas cooling, star formation and feedback, at 
small time-steps of size AT/N stsps , where AT is the time inter- 
val between outputs of the underlying cluster simulation; we adopt 
-/Vsteps = 50. Since the hydrodynamical properties of the ICM and 
galaxy positions and velocities are given by the underlying simula- 
tion, the effect of RPS is estimated for each galaxy only once per 
snapshot, considering a single instantaneous stripping event and as- 
suming that the profile of the exponential disc remains unaltered for 
R < Rsa- The mass of gas removed by RPS is transferred to the 
hot gas component of the central galaxy of its fof halo, thus becom- 
ing available to the central galaxy for cooling and star formation in 
subsequent time-steps. 

The RPS effect eats away the gas disc from the outside in. Af- 
ter the first stripping event, the remaining disc gas (if any is left) 
is assumed to form an exponential disc with the same scalelength 
as before the stripping, but sharply truncated at i? str . If in a subse- 
quent snapshot the galaxy experiences a higher RP, i? stl will then 
be smaller and the gas lost to RPS will be the gas located between 
Rsa{t n ) and R$tt(t n -i), where t n and t„_i denote the current and 
previous simulation outputs, respectively. Since in our strangula- 
tion model satellite galaxies have no infall of cooling gas, their 
gaseous discs cannot be re-built. It is also unlikely that the remain- 
ing disc gas would systematically gain angular momentum so as 
to regenerate a full exponential disc. With these considerations, if 
Rai(tn) < -Rstr(in-i), the mass of cold gas removed by RPS at 
time t„ can be shown to be 

M w {t n ) = M coM (t n ) /a(tn) ~ { a(t "- l) , (10) 

1 — Js(tn-1 ) 

where M co \i(t n ) is the mass of cold gas in the galactic disc at the 
current snapshot but before RPS acts, and f a is given by 



/.(*) 



1 + 



-Rstr(t) 
Rd 



exp (-R sa (t)/R d ) . 



(11) 



This takes into account that if all the cold gas mass is between 
R = and R = i? stl in the form of a truncated exponential disc, 
then the expression $2% for the central surface density of the cold 
gas becomes 



i(t„) 



Mcold(tn) 



2n [l-f s (t n ^)]Rl 



(12) 



Conversely, if in a subsequent snapshot the galaxy experi- 
ences a lower RP, Rm(t n ) will be larger than the radius at which 
the remaining gas is found, Rs U (t n -i)- In such a case, we set 
Rsa{t n ) = Rsa(t n —i) and no gas is removed in that time-step. 

In addition to RPS, gas could be removed from the galaxy by 
turbulent stripping, caused by the generation of Kelvin-Helmholtz 
and/or Rayleigh-Taylor instabilities at the interface between the 
cold gas and the ICM. However, we do not consider this possi- 
bility as the time-scale associated with the growth of such instabili- 
ties is generally much longer than the RPS ti me-scale, as show n by 
iMcCarthv etaU J2008h in their work (see also lFont et al]|2008l) . 

Finally, we would like to emphasize once again that RP is cal- 
culated for every satellite galaxy in the entire simulation. We do not 
restrict the calculation to galaxies residing within the most massive 
groups, and we do not turn on the calculation at a selected redshift; 
RP will act on galaxies whenever condition 10 is satisfied. 
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4 RESULTS 

4.1 Radial profiles and distributions of RP 

The procedure outlined in the previous section allows us to follow 
the evolution of the properties of the intergalactic medium within 
DM haloes, obtaining for each snapshot values of the ambient den- 
sity and relative velocity for each satellite galaxy. In this section we 
discuss the distributions and radial profiles of RP that we obtain for 
satellites in the simulated clusters, and their evolution with time. 

An interesting point is to evaluate how our RP model com- 
pares with pre vious studies which assume analytic density profiles 
for the ICM dOkamoto & Nagashim3 120031 ; lLanzoni et ai]|2005l ; 
BDL08). Hence, additionally, we estimate the RP and the proper- 
ties of galaxies assuming that the ICM gas is in hydrostatic equilib- 
rium following a NFW profile, re-scaled according to the adopted 
baryonic fraction. These profiles are calculated using the concen- 
tration parameter determined for the cluster at each redshift. We 
will refer to this as the 'analytic' method, and to our new imple- 
mentation described in Section [3~Tl as the 'gas-particles' method. 

The dependence on clustercentric distance of the ICM density, 
relative velocity (squared) of the satellites with respect to the ICM 
and the resulting RP is plotted in Fig.[2]for three different redshifts, 
where we show the results obtained for one of the G15 clusters 
(taking into account the properties of the main cluster progenitor 
when considering redshifts z > 0). In all other clusters we find a 
similar behaviour, so we choose this one as representative of the 
general trends. The median values obtained with the gas-particles 
method are shown with black points, and green shaded areas mark 
the regions enclosed by the 5th and 95th percentiles of the distribu- 
tions. In these plots we include all galaxies identified as satellites 
belonging to the cluster fof halo at each selected redshift, and the 
ICM properties are sampled at the positions of the galaxies. 

The top panels of Fig. [2] show the ICM density 
profiles obtained. As the cluster grows in mass from 
log(A/ v ir//i _1 Mq) ~ 13.7 up to 15.1 between z = 2 and 
z = 0, the median density in the cluster central region (r < 0.5r V i r ) 
increases about one order of magnitude. At the virial radius, the 
increase between z = 1 and z = is smaller. The median density 
profile decreases monotonically at the present epoch over four 
orders of magnitude, in the range < r/r vl! < 2. We compare 
these results with the ICM profiles determined for the analytic 
method (shown with dotted lines in the top panels of Fig. [2}. The 
match between the median density profile determined from the 
simulation and the NFW analytic profile is excellent at z = 0; how- 
ever, for z ^ 1, NFW profiles tend to overestimate (underestimate) 
the density for r < r\„ (r > 1.5r v jr). 

The median square velocities of the satellite galaxies relative 
to the ICM are shown in the middle panels of Fig. [2] The veloc- 
ity distribution of cluster galaxies is already established at z 1 
(at least for r < 1.5r v i r ), as can be seen from the similarities be- 
tween the radial profiles at both z = 1 and z = 0. There is a trend 
of increasing median relative velocity for decreasing radius, with a 
large scatter at all radii. At z = 2, the velocity profile is flatter and 
the median velocities are lower than for z ^ 1 . 

On any given cluster, the evolution of the RP experienced by 
satellite galaxies for z ^ 1 depends mainly on the build-up of the 
ICM density over time, as the cluster velocity profile is already 
established at z ~ 1. The bottom panels of Fig. [2] show that RP in- 
creases approximately one order of magnitude at the cluster centre 
between z = 1 and 2 = 0, consistent with the increase in ICM den- 
sity. Following the trend set by the density, the RP stays at roughly 
the same level for 1 < r/rvk < 1.5 between z = 2 and z = 1, and 



increases slightly between z = 1 and z = 0. For r > 1.5r v ir, the me- 
dian RP seems to drop between z = 2 and 2 = 1, and stays roughly 
similar for z ^ 1. This happens because, in this region, the me- 
dian density decreases over time but the mean relative velocity of 
galaxies increases, thus resulting in similar RP values. 

To determine the RP in the analytic method we also need an 
estimation of the relative velocities of galaxies. If these are not 
determined using the gas particles in the simulations, the alterna- 
tives ar e either to draw the m at random from an assumed distri- 
bution l lLanzoni et al.| [2005) or to take the velocity of the satellite 
(as given by the semi-analytic model) in the rest frame of its cen- 
tral galaxy (BDL08). If the ICM is hydrostatic this latter approach 
should be a good approximation, and this is what we use for the 
analytic method. The results of this are shown with red lines in 
the middle panels of Fig. [2] The agreement between the analytic 
and gas-particles methods is again excellent at z = 0, and also at 
z = I for r < rvir, for both the median (solid line) and the scatter 
(dashed lines). At z = 1, rest-frame velocities tend to be larger for 
r > r v i r than the velocities relative to the gas particles. 

The analytical distribution of RP that results of multiplying 
the median density given by the NFW profile by the square of the 
rest-frame velocities is shown with lines in the bottom panels of 
Fig. [2] As expected from the very good agreement between both 
approaches, both for the density and the relative velocity at z = 0, 
the RP profiles are very similar at the present epoch. The agreement 
is not so good for z ^ 1 where the shape of the RP profiles are 
similar but in the analytical model the mean values and the levels of 
the percentiles shift to higher values. According to the results of our 
self-consistent numerical approach, the analytic method appears to 
overestimate the RP over most of the range considered. 

To quantify this latter point, in Fig. [3] we plot the average dif- 
ference between the RP determined from the gas-particles and the 
analytic methods, normalized to the analytical value of RP. Results 
are shown in different colours for z = 2, 1, 0.5 and 0. The left- 
hand panel shows the result of averaging over all the simulated 
clusters, and error bars show the la cluster-to-cluster scatter. From 
this plot we can clearly see that the analytic method always over- 
estimates the value of RP for r < 1.5r v h (although at z = and 
for r > 0.4r v j, it could be said that both models agree within the 
errors). The difference between the two models increases with red- 
shift; this trend is clear within r v h , but becomes noisier outside of it 
(perhaps due to the smaller number of cluster galaxies in these re- 
gions). These trends persist if we consider separately the G14 and 
G15 clusters (middle and right-hand panels in Fig. [3] respectively). 

The distribution of RP values experienced by satellite galax- 
ies within r V ir in our simulated clusters is shown in Fig. [4] for the 
G14 clusters (top) and G15 clusters (bottom) at redshifts z = 1 
(left) and z = (right). We compare results obtained from the 
gas-particles method (histograms in black lines) and the analytic 
method (histograms with dashed lines). All the distributions can be 
well fitted by Gaussian functions, and the result of the fits we ob- 
tained for these histograms is shown in TableQ] The mean value of 
RP obtained from both the gas-particles and analytic methods are 
higher for the more massive clusters, a feature that is also present 
at higher redshifts. We mentioned before that for a given cluster the 
increase of RP for z < 1 is due to the increase in ICM density, but 
the difference between the RP values for different cluster masses is 
due to a combination of increased density and relative velocity. On 
average, the ICM density within r v - a in G14 clusters can be as high 
as 75 per cent of the density at same radii in G15 clusters. The rest 
of the difference is due to the lower relative velocity of satellites in 
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z = 2 z = 1 z = 
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Figure 2. Values of ICM density (top), galaxy velocity relative to the ICM (centre) and RP (bottom) as a function of the distance to the cluster centre, at the 
current positions of all satellite galaxies identified as members of the FOF group corresponding to one of the G15 clusters. Results are shown for three different 
redshifts: 2 = 2 (left), 2=1 (centre) and 2 = (right). Diamonds indicate median values in bins of 0.1 r/r v -„, and the shaded areas mark the regions enclosed 
by the 5th and 95th percentiles in each bin. The red dotted line in the top panels is a NFW density profile with the concentration c of the halo at that redshift. 
The solid line in the central panels indicates the median velocity of the galaxies in the rest frame of the cluster, and the dashed lines enclose the 5th and 95th 
percentiles of that distribution. In the bottom panels, the solid line indicates the median RP obtained by multiplying the density given by the NFW profile by 
the rest-frame velocity; dashed lines denote the 5th and 95th percentiles. 



the smaller clusters (within r v j r , the mean velocity in G14 clusters 
is ~25 per cent of the mean v in G15 clusters). 

The range of RP values we obtain at z = for both 
models is similar to that found by BDL08, who calculate val- 
ues of RP for clusters of masses comp arable to ours, ex- 
tracted from the |Pe Lucia & Blaizol d2007l) semi-analytic cata- 
logu e. BDL08 assume a n ICM described either by an isothermal 
or a iKomatsu & Seliakl fcOOll hereafter KS01) model, in hydro- 
statical equilibrium within a NFW halo. For the case of the isother- 
mal model, BDL08's distribution of RP is fairly skewed, with a 
sharp cut-off at higher pressures; for the KS01 model they obtain a 
much more symmetrical distribution (see their figs 1 and 2, respec- 
tively). Comparing the shapes of our local RP distributions, given 
in the right-hand panels of Fig. [4] with those of BDL08, we find 
better agreement with their results for the KS01 model. The local 
mean RP values obtained from our gas-particles method (see Ta- 
bleQ} are also in good agreement with those determined by BDL08 
with the KS01 model, who find KT 11 ' 3 dyn cm for their set of 
10 14 M Q clusters, and 10" 10 7 dyn cm -2 for their 10 15 Af© clus- 
ters. 



As already noticed from Fig. [2] we can see again in Fig.|4]that 
the distributions obtained with both methods agree very well at the 
present epoch, but not so for higher redshifts. For both sets of clus- 
ters, the mean value of RP obtained from the gas-particles method 
for galaxies within the virial radius increases approximately half 
an order of magnitude from z = 1 to the present. However, in the 
analytical calculation, the mean values at both redshifts are much 
more similar (in the case of the G14 clusters, if errors are taken 
into account then the mean values can be considered equal). This 
fact, combined with the results shown in Fig. [3] indicates that the 
RP in the clusters does not grow over time in the same way in both 
models. 

To better visualize this, we plot in Fig. [5] the mean RP for 
galaxies within r v j r as a function of redshift, for both sets of clus- 
ters. Mean RP values given by the gas-particles method are shown 
with black symbols, and the shaded areas indicate the regions 
within la of the mean. Solid lines indicate the mean values ob- 
tained from the analytical calculation, with dashed lines enclosing 
the regions within ler of the mean. For both models, the differ- 
ence between the mean RP values of both cluster sets at a given z 
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Figure 3. Average difference between the RP determined from the gas particles in the simulation and the RP calculated as the product of a NFW density profile 
and the squared velocities of the satellite galaxies in the rest frame of the clusters. Results are shown for z = 2, 1, 0.5 and 0, and for all the simulated clusters 
(left), the G14 clusters only (centre) and the G15 clusters only (right). 



Table 1. Best-fitting parameters of Gaussian functions to the RP histograms of Fig. [4] (log P mm ) and a am denote the mean and dispersion of each distribution, 
respectively. 





G14 clusters 


G15 clusters 


Model 


z 


(log P ram ) 


& ram 


(log Pram) 


(7 ram 






[ft 2 dyn 


cm -2 ] 


[ft, 2 dyn 


cm -2 ] 




1.0 


-11.64±0.08 


0.73±0.05 


-10.91±0.03 


0.76±0.02 


Gas-particles 




-11.14±0.07 


0.84±0.05 


-10.63±0.05 


0.75±0.03 


Analytic 


0.0 


-11.33±0.03 


0.87±0.02 


-10.51±0.09 


0.83±0.04 


Gas-particles 




-11.18±0.05 


0.96±0.03 


-10.47±0.07 


0.81±0.03 


Analytic 



is always about one order of magnitude, being larger in the more 
massive clusters, as expected. 

In the simulations by iRoediger & Brtiggenl bOOfj) , values of 
RP of order 10~ 12 dyn cm~ 2 were c alled weak, 10~ n dyn cm" 2 
medium and 10~ 10 dyn cm -2 strong. Roediger & Briiggen (2006) 
find that a spiral galaxy with mass ~ 2 x lO n M0 subject to strong 
RP will typically lose all of its gas within ~ 50 Myr, medium RP 
will remove approximately half of the gas within ~ 200 Myr, and 
weak RP removes relatively small amount of gas (the actual gas 
loss will depend on the structure of the gaseous, stellar and DM 
components). Fig.[5]shows that for the gas-particles model the ma- 
jority of satellite galaxies in the massive G15 clusters are experi- 
encing medium level RP already at z ~ 2, and at z < 0.5 about 20 
per cent of satellites experience strong levels of RP at any given 
time. In the case of the G14 clusters, the mean RP is between 
medium and weak for z < 1.5, and few galaxies ever experience 
strong levels of RP; only at z ~ 1 does a significant fraction of 
satellites begin to experience medium levels of RP. 

In the analytic model, at higher redshifts (z > 1), the mean RP 
for the G15 clusters reaches medium values and a significant frac- 
tion of satellite galaxies already experience medium-to-strong lev- 
els of RP at z ~ 2. A similar situation occurs for the G14 clusters, 



with the values of RP shifted to the range of weak-to-medium lev- 
els, the difference between the two models being larger for higher 
redshifts. In both models the scatter around the mean is large, span- 
ning one order of magnitude at each side of the mean. The disper- 
sion in the distributions also grows slightly with time, as can be 
better appreciated from Table [T] where we can see an increase in 
a of about 10 per cent between z = 1 and 0, for both sets of clus- 
ters. This may arise as a result of the larger number of galaxies in 
clusters at the present time (about three times more galaxies than at 

2=1). 

Although the results for both models agree at z ~ for the 
different sets of clusters, the different evolution of RP over time 
is clear. While the mean RP values obtained from gas-particles 
method increase almost exponentially with redshift, the evolution 
of the mean RP in the analytic method is flatter, especially so in 
the case of the G14 clusters. The lack of strong evolution of the 
mean RP values for G14 clusters is a consequence of the over- 
estimation of RP by the analytic model for the inner parts of the 
clusters (r < 0.5r v ; r ), even at z ~ (see middle panel of Fig. [3}. 
This is because these clusters, and their progenitors at higher red- 
shift, are galaxy group-sized systems of virial masses in the range 
~ 10 13 - 10 14 ft." 1 M Q , whose ICM distributions are not so well 
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Figure 4. Histograms of the RP experienced by the galaxies within r v j r in 
the simulated clusters at z = 1 (left) and z = (right), obtained from the 
gas-particles method (black lines) and the analytic method (dashed lines). 
Results for the G14 and G15 clusters are shown in the top and bottom pan- 
els, respectively. Error bars denote the la cluster-to-cluster scatter and for 
clarity are shown only for the gas-particles method. 



described by a NFW profile. There is also a clear but smaller dif- 
ference between the models for the case of the G15 clusters. For 
z < 0.5, both models agree fairly well; this can be due to the fact 
that these clusters are already massive systems at z ~ 2, and so 
their ICM is better described by a NFW profile (see radial profiles 
for different cluster masses in the top panels of Fig. [2}. 

In Fig. [6] we plot RP histories given by our gas-particles 
method for randomly selected satellite galaxies of type 1 (black 
dot-dashed lines) and type 2 (black solid lines). Results are shown 
for one of the G15 clusters (left) and one of the G14 clusters 
(centre). For comparison, the right-hand panel of Fig. [6] shows 
the RP evolution for galaxies within one halo with virial mass 
~ 10 13 Mq extracted from the simulation box of one of the 
G15 clusters. There are differences between the RP histories of dif- 
ferent types of satellites. For type 2 satellites, the RP typically oscil- 
lates, increasing on average; whereas for type 1 galaxies, the typical 
history is one of monotonically increasing RP. This difference in 
the RP histories is due to the way in which positions and velocities 
of galaxies are assigned in the code, being more reliable for type 1 
galaxies. As mentioned in Section l2T2l the position and velocity of 
a galaxy are assigned by tracking the most-bound DM particle of 
its host DM subhalo. In the case of type 1 galaxies this subhalo still 
exists and is affected by dynamical friction, but for type 2 galax- 
ies it has been destroyed by tidal forces and we are then tracking 
a single DM particle. The regions enclosed by the maximum and 
minimum RP values reached by all the galaxies that are present at 
z = in the selected halo are depicted by dashed lines and shaded 
areas for the analytic and gas-particles methods, respectively. Once 
again, we see differences between the models; satellites in the an- 
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Figure 5. Mean ram pressure exerted on the galaxies within r v k of the sim- 
ulated clusters as a function of redshift. Top and bottom panels show re- 
sults for the G14 and G15 clusters, respectively. Open diamonds and filled 
squares indicate the mean value of RP for the gas-particles method and the 
shaded areas denote the regions within 1 a of the mean. Solid red lines indi- 
cate the mean RP for the analytic method, and the dashed red lines enclose 
the region within lcr of the mean. 



alytical model reach higher values of RP in all three cases, and the 
difference between both methods increases with decreasing halo 
mass. 

4.2 Effects of ram pressure stripping 

After analysing the evolution of the distributions of RP values and 
the dependence on the cluster mass obtained from two different 
models, we now focus on the values given by the gas-particles 
method and explore the influence of RPS on galaxy properties. 
The new estimation of the size of galactic discs, described in Sec- 
tion [T2] affects the results of the disc instability process imple- 
mented in sag. This calls for a retuning of the free parameters of 
sag in order to retain a good fit to observed properties such as 
the bj- and it-band cluster luminosity functions and the BH-bulge 
mass relation, among others. Three parameters change with respect 
to those used by LCP08 (see the descriptions of the model param- 
eters in LCP08). Specifically, the parameters kagn and /bh, which 
control the rate of gas accretion onto the central BH, increase from 
2.5 x 10" 4 Mq yr- 1 to 10" 3 M Q yr _1 , and from 0.015 to 0.04, 
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Figure 6. RP histories of satellite galaxies in one of the G15 clusters (left), one of the G14 clusters (centre) and in one group with M v j r ~ 10 13 /i _1 Mq (right) 
which is contained in the simulation box of one of the G15 clusters. Solid black lines correspond to a selected type 2 satellite, and black dot-dashed lines to 
a type 1 satellite. Shaded areas mark the regions enclosed by the maximum and minimum RP reached by all the galaxies present in the cluster at z = as 
determined from the gas-particles method; dashed thick lines denote the maximum and minimum values of the RP calculated with the analytic method. 



respectively; the disc instability threshold edisc decreases from 1.1 
to 0.85. This recalibration is done without considering RPS. We 
then run three sets of models: one without and two with RPS, using 
the same set of parameters in all cases in order to evaluate the effect 
of RPS on galaxy properties. In the following, we will use sag to 
refer to the semi-analytic model without RPS, and sagrp and sagrp- 
a will denote the models where RP values are determined using the 
gas-particles or the analytical methods, respectively. 

As a test of the performance of the model, Fig. [7] shows 
the stripping radius of disc satellite galaxies for the sagrp model 
as a function of the peak RP that they experience. Re- 
sults for the mean stripping radius and la scatter are plot- 
ted for disc satellite galaxies in two different mass ranges: 
10 10 h _1 Af s < Afsteua < 7 x 10 10 h^Mo (dot-dashed 
lines) and 7 x 10 10 ft _1 M Q < M ste n ar < 2 x 10 11 ft _1 M 
(solid lines). We compare this with the results of sev- 
eral N-body/hydrodynamical simulations specially focused on 
the study of RPS in individual galaxies (denoted by dif- 
ferent symbols). These simulations emplo y different proce- 
dures to compute the hydrodynamics (SP H: lAbadi et al. 1999, 
Schulz & Struck! 1200 ll. Ijachvm et al.ll2007l : Euleria n: lOuilis et al. 
200( 1 iRoediger & Henslej 120051 ; sticky particles: IVollmer et al.l 
200 11) . All the simulations consider Milk y Way-class galaxies en- 
tering a dense environment, except for ISchulz & Struck fcOOll) 
who study a smalle r galaxy with M stc iiar — 6 x 10 10 h^ 1 Mq. 
Ijachvm et alj J2007h run simulations for different types of clus- 
ters; here we select the results for their standard cluster model. The 
dashed line shows the analytical estimate for R ai of Eq. $9$ for 
a Milky Wa y-like galaxy, w ith surface densities for stars and gas 
taken from Flvn n et all J2006h . There is a good agreement between 
the GG72 estimate and the different simulations (f or clarity we do 
not plot the results from Roediger & Hensler 2005, as they are very 
similar to the GG72 line). 

Even though a comparison between our method and the de- 
tailed simulations is not direct, since the initial conditions are not 
necessarily the same, the general trend seen in the detailed simu- 
lations is reproduced in our model for galaxies of similar stellar 
masses. We note that the comparison does not include those satel- 
lite galaxies from our simulations whose cold gas is completely 
stripped (i? stl = 0). There are no such galaxies for log P ram ,pcak < - 



11.5, and the majority of fully stripped galaxies (>80 per cent) are 
found at log P r am, P eak > -10.5. 

One of the galaxy properties most directly affected by RPS 
is the cold gas content. Fig. [8] shows the fraction of galaxies that 
have lost all of their cold gas (/ n o-gas) as a function of clustercentric 
distance, for G14 and G15 clusters. To construct this plot, we se- 
lect galaxies within 2r\„ of the simulated clusters and with stellar 
mass ^ 10 9 M©. For this plot all galaxies within r v ir are con- 
sidered, not only members of the main fof group; so we are includ- 
ing members of outlying groups which are infalling into the main 
cluster. The fractions of gas-depleted galaxies are determined for 
the three different semi-analytic models considered, which are rep- 
resented by different line types. Different colours identify the red- 
shifts at which results are shown. For the sag model, /nogas is very 
low; only in the innermost regions (r < 0.2r v i r ) the gas fractions 
increase from /no-gas ~ 0.2 to 0.3 in the redshift interval 1 < z < 0, 
regardless of cluster mass. These galaxies are those that have spent 
enough time within the cluster to gradually consume their cold gas 
reservoir; recall that cooling flows are suppresed in satellite galax- 
ies due to strangulation. In sag, /no-gas has an almost flat radial dis- 
tribution beyond 0.5r v h (dashed lines in Fig.[8}. 

When RPS is acting, almost all galaxies in the innermost re- 
gions of clusters lose their cold gas by z = 0, for both sets of clus- 
ters. The fractions of gas-depleted galaxies in the sagrp model de- 
crease for larger distances to the cluster centre, but they remain 
significant at r ~ 2r v h, being larger for G15 clusters (/no-gas — 0.5) 
than for G14 ones (/no-gas — 0.3). The lower values of / no -gas in the 
outskirts of G14 clusters in the models that include RPS are also 
consistent with the weak-to-medium values that RP takes in these 
clusters. This monotonic decrease of /no-gas with clustercentric dis- 
tance is a direct consequence of the behaviour of the radial profile 
followed by RP values (see bottom panels of Fig.[2jl. 

The evolution of the cold gas content of galaxies in the 
sag model is very similar for clusters of different masses; however, 
this is not the case for galaxies in models which include RPS. The 
difference of over half an order of magnitude between the mean 
RP values for G14 and G15 clusters at different redshifts, which 
is clearly visible in Fig. [5] drives a rather different evolution of 
the fraction of galaxies devoid of cold gas. In the sagrp model, for 
G15 clusters the situation at z = regarding the cold gas content of 
galaxies has been practically established at z = 0.5, and it does not 
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evolve much since 2 = 1. Although at this latter redshift the frac- 
tions are m 10-20 per cent lower than at z = at all radii, they are 
quite high (/no-gas > 0.5) even in the cluster outskirts. Conversely, 
for the case of the G14 clusters we see a stronger evolution in the 
fractions of gas-depleted galaxies in the same redshift interval. This 
reflects a more gradual removal of the cold gas of satellites as a re- 
sult of the lower values of RP in these clusters. In the outskirts of 
the G14 clusters / no -gas ~ 0.3, lower than for G15 clusters but still 
much larger than in the sag model. 

The results described above clearly show a dependence of the 
evolution of the fraction of gas-depleted galaxies on cluster mass. 
This trend is different for the analytic estimation of RP used in the 
sagrp-a model, shown by dotted lines in Fig. [8] Both for the G14 
and G15 clusters, the final values of /no-gas are reached at higher 
redshift than for the sagrp model. For instance, in the case of the 
G14 clusters, the value of /no-gas in the sagrp-a model at z = 1 is 
already as high as the value reached for the sagrp model at z = 0.5. 
This behaviour is the result of the much higher values of RP at 
z > in the sagrp-a model (see Fig.|5j. 

In the models that include RPS, the higher fractions of gas- 
depleted galaxies found at z = in the cluster outskirts, for both 
sets of clusters, indicate that during cluster assembly RPS has an 
important effect on the galaxies contained in the smaller subgroups 
that are being accrete d by the cluster; this is wh at has been called 
'pre-processing' (e.g. Fuiita 20041; lMihosll2004l) . This idea is also 
inferred from observational evidence of compact group of galax- 
ies falling into the massive galaxy clusters Abell 1689 (z ~ 0.18), 
Abell 2667 Q ~ 0.23) and Abell 1367 (z ~ 0.02) JCortese et al] 
I2006L l2007h ; galaxies within these groups show extended tails of 
ionized gas, that can be interpreted as the signature of RPS. We 
then find that RPS is a process that significantly affects the cold 
gas content of galaxies and the rate at which it is removed, which 
depends on cluster mass. Consequently, the RPS experienced by 
satellite galaxies might enhance differences in the star formation 
history and the colour evolution on haloes of different masses; this 
will be explored in detail in a forthcoming paper. 



5 CONCLUSIONS 

We develop a new method to describe the effects of RPS on galaxy 
properties which works within the sag semi-analytic model of 
galaxy formation and evolution, in combination with cosmologi- 
cal hydrodynamical iV-body/SPH simulations of galaxy clusters. 
RPS is implemented in sag adopting the GG72 criterion. The novel 
feature of our implementation is that the kinematical and fhermo- 
dynamical properties of the hot gas responsible for RP are ob- 
tained from the gas particles of the SPH simulations (gas-particles 
method). This results in a more self-consistent estimation of the RP 
experienced by satellite galaxies. 

We compare our results with those obtained from an analytic 
estimation of RP, which considers a NFW density profile for the hot 
gas contained within DM haloes, re-scaled according to the adopted 
baryonic fraction (analytic method). We analyse the dependence on 
clustercentric distance and redshift of the RP values given by both 
methods, and evaluate the influence of the environment on the be- 
haviour of these distributions. The RPS method discussed in this 
work can be adapted to work with any kind of cosmological simu- 
lation that includes gas physics, regardless of the particular numer- 
ical scheme used for the hydrodynamical calculations. However, 
the specific results from its application might depend on the details 
of the numerical implementation. For example, the SPH method 




K111 1 1 1 \, 1 

-14 -13 -12 -11 -10 -9 
log P ram , P eak [h 2 dyn cm 2 ] 



Figure 7. Stripping radius i? str as a function of the peak RP experi- 
enced by disc satellite galaxies (Afbulge/A^stellai ^ 0.95). Thick dot-dashed 
and solid lines denote the mean Rsa for satellites with 10 10 h~ 1 MQ < 
Miliar < 7 x W 10 h- 1 M Q and with 7 x 10 10 h -1 M Q < 
A^stellar < 2 X 10 11 h~ 1 Mq, respectively, and thinner lines enclose the 
regions within la of the mean values. Points show the results of differ- 
ent N-body/hydrodynamical simulations of RPS in single galaxies. The 
dashed line shows the analytical estimate of equation |9) for a Milky Way- 
class g alaxy, with surface densities of stars and gas taken from Flynn et al. 
120061) . 

utilises an artificial viscosity term to properly capture hydrody- 
namic shoc ks, and this viscosity can artificially su ppress turbulence 
in the ICM toolag et alj2005l ; lAgertz et al.l2007h . The simulations 
we have used include the standard SPH artificial viscosity. The use 
of a lower-viscosity formulation of SPH, or a grid-based method, 
may result in an increased level of ICM turbulence. This could re- 
sult in the distribution of the velocities of galaxies relative to the 
gas deviating further from the velocities calculated by assuming a 
hydrostatic ICM. 

We have selected from the simulations considered two sets of 
cluster-sized haloes, with masses Mvir ~ 10 14 /i _1 M© (G14 clus- 
ters) and M yil ~ 10 15 ft _1 Af (G15 clusters). The RP values ob- 
tained by the gas-particles and analytic methods are used by the 
RPS process implemented in the semi-analytic model, resulting in 
the models sagrp and sagrp-a, respectively. These models provide a 
galaxy population affected by RPS; their results are compared with 
those obtained from the standard sag, which does not take RPS into 
account, in order to evaluate the effect of RPS on the cold gas con- 
tent. We summarize our main results. 

(i) The RP estimated from the gas-particles method increases 
approximately one order of magnitude at the cluster centre be- 
tween z = 1 and 0, consistent with the increase in ICM density, 
since the cluster velocity profile is already established at z = 1. 
Median RP values do not evolve much in the outskirts of the clus- 
ter (r/r v i r > 1), which is a consequence of the combined behaviour 
of the median ICM density and the relative velocities, whose radial 
distributions decrease and increase with time, respectively. 

(ii) The radial distributions of median ICM density, velocity rel- 
ative to the ICM and RP in the gas-particles method follow smooth 
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Figure 8. Fractions of galaxies that have lost all of their cold gas as a function of distance to the cluster centre, for G14 clusters (left) and G15 clusters (right). 
Results are shown for three different redshifts: 2=1, 0.5 and 0, and for three models: sag (dashed lines), sagrp (solid lines), which takes into account RP 
values estimated from the gas-particles method, and SAGRP-A (dotted lines), where the RP is calculated with the analytic method. All galaxies within 2r v - a of 
the simulated clusters and with stellar mass greater than 10 9 /i" 1 Mq are included in this plot. Galaxies are binned in 0,4r/r v ir, and error bars are only shown 
on the SAGRP model for clarity. 



profiles at z = 0, and these distributions are very well traced by 
analytical estimations based on NFW profiles for the ICM density. 
However, the agreement is not so good for z ^ 1, where the shape 
of the RP profiles are similar but the mean values and levels of 
the percentiles shift to higher values in the analytical model. We 
find that the analytical calculation systematically overestimates the 
RP when compared to our self-consistent numerical approach. This 
overestimation grows with increasing redshift. 

(iii) The distribution of RP values experienced by satellite 
galaxies within r v ir, estimated by both the gas-particles and the an- 
alytic methods, can be well fitted by Gaussian functions. The mean 
values of RP obtained from both methods are higher for the more 
massive clusters, a feature present for all redshifts. However, the 
difference between the two models becomes larger as one goes to 
higher redshifts, and with decreasing halo mass. This happens be- 
cause the ICM density profiles in less massive clusters, and their 
progenitors at higher redshift, are galaxy-group sized systems of 
virial mass ~ 10 13 — 10 14 ft" 1 M , whose ICM distributions are 
not so well described by the analytical profile adopted. In the ana- 
lytical model the evolution of the mean RP with redshift is milder, 
especially in the case of the G14 clusters. 

(iv) For the gas-particles method, the RP values at z = are 
weak (4.68xlCT 12 h 2 dyn cm" 2 ) and medium (3.09x 1CT 11 h 2 
dyn cm -2 ) for G14 and G15 clusters, respectively. The majority of 
satellite galaxies in massive G15 clusters are already experiencing 
medium-level RP at z ~ 2, and at z < 0.5 about 20 per cent of satel- 
lites experience strong levels of RP (of order 10 -10 dyn cm -2 ). In 
the case of the G14 clusters, the mean RP is between medium and 
weak for z < 1.5, and few galaxies ever experience strong levels of 
RP; only at z ~ 1 does a significant fraction of satellites begin to 
experience medium levels of RP. 



(v) Based on our self-consistent approach, we find that RPS has 
a strong effect on the cold gas content of galaxies in both sets of 
clusters. At z = most galaxies (> 70 per cent) within r v j r are 
completely depleted of their cold star-forming gas. This is a strong 
difference with the model without RPS, where most galaxies man- 
age to retain some cold gas: at the present epoch, only in the clus- 
ter cores (r < O.Srvir) the fractions of gas-depleted galaxies reach 
~ 40 per cent. Observations of gas fractions as a function of clus- 
tercentric distance could provide strong constraints for the models. 

(vi) The rate at which the cold gas is stripped from satellite 
galaxies depends on the virial mass of their host clusters. In our 
sagrp model, the fractions of gas-depleted galaxies for G 14 clusters 
increase appreciably between z = 1 and 0, whereas for G15 clusters 
the fractions at the present epoch are mostly established already at 

2=1. 

The general picture that emerges from the main results sum- 
marized above are that the RPS effect depends on halo virial mass 
and redshift, being more important in more massive haloes. Less 
massive galaxies within larger haloes are the most affected, so this 
could be the mechanism responsible fo r the transformation of dlrr 
galaxies into dSph in galaxy clusters jBoselli et alj|2008h . In the 
more massive clusters the gas removal is extremely effective. RPS 
could contribute to the pre-processing of galaxies in smaller groups, 
before they fall into larger, cluster-sized systems. 
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